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I.  SUMMARY 


This  is  •  User's  Manual  for  the  computer  program  UCIN-CABLE  IZ.  The 
program  is  designed  to  study  the  three  dimensional  dynamics  of  submerged 
and  partially  submerged  towing  cables.  The  cables  themselves  may  have 
multiple  branches.  Figure  1.  shows  a  schematic  illustration  of  a  typical 


cable  system.  The  system  can  of  course  be  much  simpler  than  that  shown  in 
Figure  1.  For  example*  It  could  consist  of  a  single  cable  segment  as  with  a 
buoy  or  balloon  mooring  system  as  depicted  in  Figure  2. 


Buoy  or 
Balloon 


Figure  2.  A  Simple  Mooring  System 
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The  manual  provides  Instruction  for  using  CABLE  to  study  sutri  cable 
systems.  It  also  provides  sample  Input  and  output  data.  The  computer 


program  Itself  Is  written  In  FORTRAN.  Hence*  the  Input  requirements  are 
In  the  FORTRAN  format. 


(  II.  CAPABILITIES  OF  CARLE 

CABLE  1$  designed  to  perform  the  following  dynamic  analyses: 

Given:  a)  the  physical  dota  of  the  cable  (weight*  diameter,  length)  and 
the  physical  data  of  the  towed  body;  b)  the  configuration  of  the  cable 
branches;  c)  the  fluid  properties;  ti)  the  motion  of  the  towed  end  of  the 
cable;  and  e)  the  Initial  configuration  of  the  cable  system;  then  the 
program  determines  the  kinematics  (position,  velocity,  and  acceleration) 
at  the  various  points  of  the  cable  system  as  well  as  the  cable  tension.  ^ 

i  /' 
/ 

CABLE  provides  the  user  with  the  following  options: 

1)  Either  English  cr  metric  units  may  be  used. 

2)  A  two-  or  three-dimensional  analysis  may  be  selected. 

3)  The  cable  system  may  be  Immersed  In  two  fluid  media  - 
normally  air  and  water. 

4)  The  fluid  media  may  be  given  a  uniform  "strear"  velocity. 

5)  Fluid  forces  may  be  applied  to  the  cable  system.  These 
forces  Include: 

a)  Normal  drag  forces 

b)  Tangential  drag  forces 

c)  Added  mass  forces  due  to  cable  acceleration 

d)  Bouyancy  forces 

Theso  fluid  forces  may  be  partially  or  totally  neglected. 

6)  The  towed  end  of  the  cable  may  have  arbitrary  motion  relative 
to  a  ship  frame.  The  ship  frame  Itself  may  have  motion  In 

a  plane. 

7)  Gravity  forces  may  be  Included  or  neglected. 
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II.  Capabilities  of  CABLE 


In  addition  to  the  above  or' Iona,  the  uaer  may  arbitrarily  select 
the  numerical  Integration  parameters.  Certain  output  options  are  also 
available. 


III.  THEORETICAL  BASIS  OF  CARLE 


The  cable  system  Is  modelled  by  a  series  of  cylindrical  links. 
These  links  are  connected  In  a  chain  to  simulate  the  cable.  Figure  3. 
depicts  such  a  modelling.  The  number  of  links  and  hence  their  length. 
Is  a  user  option. 


Figure  3.  A  Finite-Segment  Model  of  a  Cable. 

The  Rodel  of  Figure  3.  Is  a  "finite-segment"  model  of  the  cable. 

As  such  It  forms  a  "general-chain"  system  or  an  "open-tree"  system. 
References  Cl— 33*  provide  an  analytical  basis  for  studying  these  systems. 
The  governing  differential  equations  are  obtained  using  Lagrange's  form 
of  d'Alembert's  principle  and  Kane's  equations  as  exposlted  by  Kane  and 
oilers  E4-83.  This  procedure  has  distinct  advantages  over  Newton's  laws 
and  Lagrange's  equations  for  analyzing  multibody  systems.  Specifically, 
Kvne's  method  provides  for  the  automatic  elimination  of  non-working  Internal 
constraint  forces  while  avoiding  the  tedious  differentiation  of  scalar 
energy  functions. 

In  the  formulation  of  CABLE  the  velocities  and  accelerations  are 
computed  through  vector  derivatives  which  may  be  evaluated  through  vector 
cross  products.  This  leads  to  algorithms  which  are  readily  converted 
Into  computer  subroutines. 


*  Numbers  In  brackets  refer  to  References  at  the  end  of  the  manual. 
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III.  Thtgrttlcjl  Basts  of  CABLE 


Tht  fluid  fore**  art  computed  using  the  procedures  outlined  In 
References  119,103*  These  Include  Inertia  ("added  ness")*  dreg  (normal 
and  tangential)*  and  hydrostatic  (buoyancy)  forces.  They  ire  represented 
o>'  the  cable  links  by  forces  passing  through  the  centers  of  the  links 
together  with  couples. 

Finally*  the  governing  equations  are  solved  using  RK6S  C113*  a  fourth- 
order  Runge-Kutta  numerical  Integrator.  CABLE  Is  written*  however*  so 
that  other  numerical  Integrators  may  also  be  used. 


IV.  DEFINITIONS  OF  TERMS 


This  part  of  the  Manual  defines  the  terminology  used  In  CARLE. 
These  definitions  are  useful  In  understanding  the  Input  requl rements. 
Under 1  .  ng  Is  used  to  further  Identify  terms  of  Interest. 

1.  Cable  Model 

Figure  4.  shows  a  finite-segment  cable  model  of  a  cable  system. 


Figure  4.  A  Finite-Segment  Model  of  a  Cable  System. 

It  consists  of  a  series  of  cylindrical  cable  links  connected  In  "tree-like" 
fashion  such  that  no  closed  loops  are  formed. 

2.  Cable  Link  System 

A  cable  link  system  Is  a  set  of  connected  cable  links  together  with 
the  towed  bodies  as  shown  In  Figure  4.  Thus*  a  cable  link  system  Is  a 
multibody  system  as  defined  In  References  C2,33. 

3.  Cable  Link 


A  cable  link  Is  an  Individual  member  of  a  cable  link  system  (aside 
from  towed  bodies).  A  cable  link  Is  a  rigid  cylindrical  rod.  It  Is 
connected  to  adjacent  links  and  the  towed  bodies  by  either  spherical  or 
hinge  joints. 


IV.  Definitions  of  Ttnii 


4.  Reference  Links  and  Reference  Point! 


The  uppermost  1  Ink  of  a  cable  link  system  Is  called  the  reference  link 
*s  shown  In  Figure  5.  The  upper  end  of  the  reference  link  is  called  the 
system  reference  point  Q.  In  a  typical  dynamic  configuration  of  the  system, 

Q  Is  given  a  prescribed  motion  relative  to  the  mean  ship  frame  (See  7.  below.). 

Also,  each  of  the  links  of  the  cable  link  systems  has  Its  own  reference 
point.  For  the  reference  link.  It  Is  Q  the  system  reference  point.  For  any 
other  link,  say  Lfc  ,  It  Is  the  joint  location  at  the  upper  e>nd  of  the  link 


Reference  Link 


Figure  5.  Reference  Link,  Reference  Point,  and  Towed  Bodies 


5.  Towed  Bodies 


Currently  In  CABLE  there  Is  provision  for  towed  spherical  bodies  at 
the  ends  of  the  branches  of  the  cable  system.  These  bodies  are  referred 
to  as  the  towed  spheres.  The  towed  spheres  are  attached  to  the  cable  links  by 
by  either  a  spheri . ^1  or  hlnoe  joint  at  the  center  of  the  sphere. 
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Link  Connection  A rra> 


The  conflagration  of  the  cable  link  system  and  the  arrangement  of 
the  branches  can  be  described  by  a  link  connection  array*  To  develop  this 
array,  let  the  links  of  the  system  be  nuabered  as  follows:  Let  the 
reference  link  be  link  1,  called  L^«  Next,  numbar  the  remaining  links  and 
towed  spheres  In  ascending  progression  away  from  Lj  through  the  branches 
of  the  system.  Finally,  let  the  mean  ship  frame  be  body  0,  called  LQ. 
Figure  6.  Illustrates  this  numbering  procedure  for  the  cable  system  of 
Figures  4.  and  5. 


Figure  6.  Numbering  of  the  Cable  Link  System  of  Figures  4.  and  5 


XV.  Definitions  o  f  Terms 


If  the  links  and  towed  spheres  arc  numbered  In  this  manner*  each 
link  and  sphere*  except  the  reference  link*  Is  connected  to  one  and  only 
one  adjacent  lower  numbered  link.  (Note*  however*  that  a  link*  for 
example,  may  be  connected  to  more  than  one  adjacent  higher  numbered 
link.) 

\ 

The  array  listing  the  lower  numbered  links  Is  called  the  link 
connection  array.  For  the  system  shown  In  Figure  6.  this  array  Is: 

I 

j  Link/Sphere:  1  2  3  4  5  6  7  8  9  10  11  12  13  14 

Array:  01  2  345  2  78  9  8  11  12  13 

j 

Note  that  the  configuration  of  the  system  can  be  constructed  once  the  link 
connection  array  Is  known.  That  Is*  there  Is  an  equivalence  between  the 
i  cable  system  configuration  and  the  connection  array. 

7.  Kean  Ship  Frame*  Inertia  Frame,  and  Coordlncte  Systems 

In  CABLE  the  reference  point  9  at  the  upper  end  of  the  reference  link 
may  be  given  an  arbitrary  prescribed  motion  relative  to  a  wean  ship  frame 
which  Itself  may  be  given  a  prescribed  motion  relative  to  an  Inertial  refer¬ 
ence  frame.  This  Is  depicted  In  Figure  7.  The  mean  ship  frame  Is  assumed 


J 


IV.  Definitions  of  Term 


to  be  rigidly  attached  to  a  towing  vessel  such  as  a  ship  or  helicopter. 

Coordinate  system  are  associated  with  both  the  Inertia  fram  and  the 
man  ship  fram  as  shown  In  Figure  7.  In  the  man  ship  fram  the  Xg  dli-ection 
Is  forward"  Zg  Is  vertically  down  and  Yg  Is  In  the  starboard  direction. 
Although  the  mtlon  of  Q  relative  to  the  man  ship  fram  my  be  arbitrarily 
specified*  the  notion  of  the  man  ship  fram  relative  to  the  Inertia  fram 
Is  restricted  to  notion  along  Xg  and  rotation  (steering)  about  Zg. 

Finally,  CABLE  has  the  option  of  locating  the  man  ship  fram  at  at 
elevation  h  (which  could  be  negative)  above  a  fluid  Interface  or  water 
level. 

8.  Acceleration  Profiles 


The  specification  of  the  notion  of  reference  point  ft  relative  to  the 
man  ship  fram  and  the  notion  of  the  man  ship  fram  relative  to  the 
Inertia  fram  (See  Figure  7.)  nay  be  aeconpllshed  either  by  using  an 
acceleration  profile,  described  here,  or  by  using  a  precoded  function. 
described  In  the  following  section. 

An  acceleration  profile  Is  tlnply  a  set  of  data  representing  the 
coordinates  of  selected  points  of  an  aceele rat Ion-tine  curve.  Such 
coordinates  nay  be  obtained  froe  the  graph  of  the  acceleration  function. 

CABLE  has  the  caoablllty  of  accepting  as  a any  as  25  data  points  fron 

an  acceleration  curve.  A  piecewise  linear  approxlnation  1$  then  nade  of 

the  acceleration  function.  For  exanple,  consider  the  acceleration  function 

and  Its  approxlnation  shown  In  Figure  8.  The  acceleration,  velocity,  and 

^  h 

dlsplacenent  during  the  1  tine  Interval  are  then: 
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IV.  Definition!  of  Tersr 
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Figure  8.  Acceleration  Profile  Approximation 

a  *  a .  +  (  *~1 — *1  )(t  -  t  )  (1) 

1  t1+1  "  *1  1 

(2) 

t,)3^  «> 

where  a^»  v^,  d^,  and  t^  are  the  acceleration,  velocity,  displacement, 
and  time  at  the  beginning  of  the  1*^  Interval.  Thus,  the  entire  kinematic 
profile  (displacement,  velocity,  and  acceleration)  is  known  when  the  a^ 
are  given  and  when  v^  and  d^,  the  Initial  velocity  and  displacement 
(at  time  t . ’ ,  are  given. 

9.  Precoded  Functions 

If  the  motion  of  reference  point  fl  relative  to  the  mean  ship  frame 
or  If  the  motion  of  the  mean  ship  frame  relative  to  the  inertia  frame  is 


d,  +  v,(t  -  t.)  +  a.(t  -  t,)V2  +  (  ' 
1111  1  t.+1  -  t1 
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IV.  Definitions  of  Term* 


relatively  simple.  It  may  be  convenient  to  describe  the  motion  using  a 
precoded  function.  For  example,  these  functions  may  be  used  to  describe 
constant  or  decaying  sinusoidal  speed  of  the  mean  ship  frame.  Sped fl ally, 
CABLE  provides  the  user  with  the  following  precoded  function: 

f  *  f q  +  Ae***  cosCpt  +  4)  (4) 

where  fQ,  A,  b,  p,  and  4  are  user  supplied  constants.  Hence,  for  a  constant 
ship  speed,  the  user  would  simply  provide  a  value  for  fQ  with  the  other 
values  being  zero  (which  are  the  default  values).  Specific  Input  format  Is 
discussed  In  the  following  part  of  the  Manual. 

10.  Fluid  Velocities 


CABLE  also  has  the  provision  of  allowing  each  of  the  fluids  to  have 
a  uniform  (constant  In  time  and  space)  "stream'1  velocity.  Specific  Input 
format  is  described  In  the  following  part  of  the  Manual. 

11.  Link  Orientation  A  qles 

/ 

The  relative  orientation  of  two  adjacent  links  can  be  defined  In  terms 
of  their  relative  orientation  angles  as  follows:  Consider  a  typical  cable 
link  L^  as  shown  In  Figure  9.  Let  X^,  Y^,  and  represent  a  coordinate 


Figure  9.  Typical  Link  L^  and  Coordinate  Axes. 
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IV.  Definitions  of  Terms 


frane  attached  to  L^  such  that  the  axis  of  Is  along  X^  as  shown.  Next-, 
consider  link  and  Its  adjacent  lower  numbered  link  as  shown  In 
Figure  10.  Let  the  axes  system  X^,  Yj,  and  of  be  Initially  aligned 


fLk  ^LJ 


Figure  10.  Two  Typical  Adjoining  Links  and  Aligned  Axes  Systems. 

with  the  axes  system  of  L^  as  shown.  Now*  Lk  can  be  brought  Into  a  general 
orientation  relative  to  L^  by  three  successive  rotations  through  angles 
Yk*  ak*  *nd  Sk  •bout  axe*  Z|c*  \*  and  Y|(*  The**  angles  are  taken  as 
positive  when  the  rotation  Is  "right-handed"  or  dextral  relative  to  the 
axis.  They  are  thus  sometimes  called  "dextral"  orientation  angles.  For 
example,  a  positive  Yk  rotation  Is  shown  In  Figure  11. 
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TV.  Definitions  of  Terms 


Not*  that  th*  angles  y^,  a^»  and  0k  art  defined  In  teras  of  rotation  of 
about  the  Z^,  X^,  and  Yk  axes  (as  opposed  to  rotations  about  Zj,  ,  and 
Yj).  This  aeans  that  If  yk  **  0,  ak  1o  the  rotation  of  about  the  "rotated" 

axis  as  shown  In  Figure  1?.  6^  Is  defined  similarly  as  a  rotation  about 

th*  rotated  Yfc  axis. 


Figure  12.  Positive  Rotation 

The  rewson  for  th*  ordering  yk,  a^,  0k  (as  opposed  to  say  a^,  0k,  yk> 
is  to  have  the  "middle"  rotation  (that  Is,  ak>  to  be  about  the  cable  link 
axis.  This  alddle  rotation  Is  then  a  aeasure  of  the  cable  "twist."  Assualng 
that  the  twist  between  adjacent  links  Is  relatively  small  (specifically, 
less  than  90°),  this  ordering  (y^,  ok»  0k>  avoids  orientation  singularities 
as  discussed  In  Reference  C123. 


Initial  Cable  Configuration 


The  Initial  configuration  of  the  cable  system  before  the  equations  of 
aotlon  are  Integrated,  depends  i,pon  the  Initial  values  of  the  relative 
orientation  angles  of  th*  links  and  their  derivatives.  For  example.  If 
the  Initial  values  of  all  the  orientation  angles  are  zero,  the  Initial 
configuration  of  a  single  branch  cable  system  would  be  as  shown  In  Figure  13. 


IV.  Definitions  of  Terns 


13.  Two  and  Thro  D1 Mansions!  Analyses 


CABLE  ha*  baan  developed  to  provlda  a  three-dlnenslonal  analysis  of 
cable  dynanlcs.  However,  thara  ara  many  eaaaa  of  Interest  whara  tha 
cabla  aovas  nearly  In  a  plana.  To  acconodate  tha*#  casts,  CARLE  has  a 
two-dlnenslonal  analysis  option. 

14.  Spherical  and  Hlnoa  Connacting  Joints 

If  tha  cabla  has  thraa  dlnenslona!  Motion,  tha  cabla  links  art  connactad 
to  aach  othar  by  spharlcal  (that  Is,  "ball-and-socket")  joints.  Whan  the 
notion  Is  restricted  to  a  plana  as  In  tha  two-dl Mansions!  option,  the  cable 
links  are  connactad  to  aach  othar  by  pin  (that  Is,  "hinge")  joints.  In 
this  latter  case,  only  one  of  tha  thraa  orientation  angles  (0|<)  Is 
needed  to  define  tha  orientation  of  a  typical  link. 

15.  Weight  Density,  Link  Length  and  Link  Planetar 

The  physical  paraneters  needed  by  CABLE  to  describe  the  cable  links 
are  tha  weight  density,  link  length,  and  link  dlaneter.  The  weight  density 
Is  defined  as  tha  weight  par  unit  length. 

Tha  cable  links  nay  have  different  lengths,  dlsMeters,  and  weights. 

That  Is,  tha  cabla  need  not  be  unlfom  In  Its  physical  properties. 

What,  the  above  data  ara  given,  CABLE  autonntlcally  calculates  the  link 
nass  canter  locations  and  the  link  Inertia  dyadic*. 

16.  Weight  and  Planetar  of  tha  Towed  Spheres 

In  addition  to  the  link  physical  paraneters,  CABLE  also  requires  the 
physical  paraneters  of  the  towed  spheres.  These  are  slnply  the  sphere 
waists  and  dlaneters. 


IV.  fteflnjtjona  of  Terms 


17.  Fluid  Forces 


CABLE  hat  tha  options  of  Including  up  to  four  kinds  of  fluid  forcas 
exart  ad  on  tha  cable  links?  1)  drag  forcas  normal  to  tha  link?  2)  drag 
forcas  parallel  (tangent)  to  tha  link)  Ji  "addad-aass"  forcas)  and 
4)  buoyancy  forcas.  Rafaranca  C93  describes  these  forcas  In  detail. 

Tha  drag  forces  and  tha  added  aass  forcas  may  be  defined  as  follows: 
Consider  a  typical  cable  link  as  shown  In  Figure  14.  Let  6^  be  tha.  canter 


Figure  14.  Fluid  Force  at  a  Typical  Point  of  a  Cabla  Link, 
of  tha  link  and  let  P  be  a  point  a  distance  x  frow  6^  as  shown.  Tha  added 


aass  and  drag  forcas  at  P  aay  than  be  represented  as: 

t  *  a!n  +  bI*nW  *  cI*tI*t  (5) 

where  yN  Is  tha  normal  component  of  tha  fluid  velocity  relative  to  tha 
cable  link  at  Pt  and  the  coefficients  A,  B,  and  C  are: 

A  «  Cjjp(ir/4)d^  (6) 

B  ■  C^p(d/2)  (7) 

C  *  Cp0(d/2)  (8) 


IV*  Definitions  of  Tera* 


where  pis  the  fluid  mss  density  and  d  Is  tha  disaster  of  the  cable  link. 
Cfl*  C^t  and  Cy  are  coefficients  dependent  upon  the  Reynolds  nuaber  of  the 
fluid  floe  past  the  cable  link.  These  coefficients  are  usually  date rained 
experl  pent  ally.  In  CABLE  the  following  values  are  used  C9»  133  s 


Si" 


1.0 


(9) 


0.0  for  Re^  0. 1 


0.45  ♦  5.93/ (Re,,)0* 33  for  0.1  <  ReN  <  400.0 
1.27  for  400.  <  ReN  <  105 


(10) 


0.3  for  Re^  >  103 

0.0  .  for  Rey  *  0.1 
1.88/(ReT)0,74  for  0.1  <  RSy  <  100.55 


0.062  for  ReT  >  100.55 


where  the  Reynolds  nuabers  Re^  and  ROj  are  defined  as 


(11) 


R*H  *  pdI^N*/w 


(12) 


and 


Ray  »  pd|</T|/y  (13) 

where  u  Is  the  fluid  viscosity. 

In  Equation  (5)  the  first  tera  Is  called  the  "added-aass"  force*  the 
second  tera  Is  called  the  "nonaal  drag"  force*  and  the  third  tera  Is 
called  the  "tangential  drag"  force.  In  CABLE  the  systea  of  added  aass 
and  drag  forces  acting  at  all  the  points  of  the  cable  link  are  replaced 
by  a  single  force  passing  through  6^  together  with  a  couple.  For  additional 
Inforaatlon*  see  References  C9*133. 
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IV.  Definitions  of  Terms 


The  buoyancy  forces  art  due  to  tha  hydrostatic  praatura  forcaa  exerted 
on  tha  cabla  link.  Thaaa  forcaa  ara  normal  to  tha  aurfaeaa  of  tha  link 
which  la  axpoaad  to  tha  fluid.  Slnca  tha  arda  of  tha  cabla  link  ara 
not  axpoaad  to  tha  fluid  tha  raaultant  hydrostatic  force  on  tha  link  la 
normal  to  tha  link  axis.  Hence,  the  system  of  hydrostatic  forces  may  be 
replaced  by  a  single  force  normal  to  tha  link  axis  given  by: 

*  -pgv^n  x  (k  x  p)  C14> 

where  o  Is  the  fluid  mass  density,  g  Is  the  gravity  constant,  Is  the 

cable  link  volume,  n  Is  a  unit  vector  along  tha  link  axis,  and  k  Is  a  vertical 

unit  vector  as  shown  In  Figure  14.  (For  additional  Information,  sea  Reference  C93.) 

Tha  fluid  forces  on  tha  towed  spheres  ara  obtained  similarly.  See 
Reference  C93. 

18.  Height  Forces 


The  weight  forces  may  be  represented  on  each  cable  link  by  a  single 
vertical  (downward)  force  passing  through  tha  mass  center  6^  of  L^.  If 
Is  the  weight  density  per  unit  length  of  l*.  Uk  nay  be  expressed  as: 


“k 


rk4k! 


(15) 


where  Is  the  length  of  and  where  k  Is  the  vertical  unit  vector  shown 
in  Figure  14. 

Similarly,  the  weight  forces  on  the  towed  spheres  may  be  represented 
by  single  vertical  forces  passing  through  the  center  of  the  spheres  with 
magnitude  equal  to  the  sphere  weight. 


19.  Units 

CABLE  allows  the  user  to  select  either  English  or  metric  units  for 
the  Input  and  output  data.  If  English  units  are  selected  the  units  are 


-19- 


IV.  btfinitions  of  Ttrw 


slugs.  fast.  and  tacondt  for  aatt,  langth,  and  tint.  Howavar,  inchat  ara 
us ad  for  tha  link  diaaatar  aid  towad  tphara  diaaatar.  If  aatric  units 
ara  salactad  tha  units  ara  kilooraas.  oat art,  and  saconds  with  cant i attars 
ustd  for  link  diaaatar  and  towad  sphara  diaaatar. 

20.  L aba  Is 


CABLE  allows  tha  uaar  to  arbitrarily  labal  or  naaa  thv  links,  thalr 
rafaranca  points,  and  tha  towad  spharas. 
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y.  INPUT  DATA 


This  pe*t  of  the  Manual  describes  the  specific  Input  data  require¬ 
ments  of  CABLE.  The  data  Itself  la  dascHbad  In  terms  of  "line  Images" 
such  at  would  also  appaar  on  a  computer  card. 

Aa  noted  earlier,  the  algorithms  and  coda  of  CABLE  are  wHttan  In 
FORTRAN.  Hence*  the  Input  data  la  to  be  submitted  In  tha  FORTRAN  formal. 

Integer  data  la  to  be  enured  In  IS  format  and  real  data  In  either  F10.4 
or  E10.4  unlec*  otherwise  specified.  The  order  of  the  Input  data  Is  In  the  same 
order  aa  It  Is  described  below.  The  units  of  the  Input  data  are  either  English 
or  metric  according  to  the  option  selected. 

1.  Heading  or  Title 

The  first  line  of  Input  to  CABLE  contains  a  title  which  may  be  used  to 
Identify  the  particular  data  entered.  It  may  contain  up  to  80  characters, 
and  it  will  appear  at  the  top  of  each  computer  output  data  page. 

2.  Options 

CABLE  has  a  number  of  user  options  to  define  the  desired  analysis. 

7h«s*  options  are  selected  by  specifying  Integer  values  (usually  0,1,2,3,4, 
or  5)  a*  follows: 

1)  Units  tSee  Section  19.,  Part  IV):  On  the  first  line  of  the 
data  enter  a  1  If  English  units  are  desired.  Enter  a  2  If 
metric  units  are  desired. 

2?  Two  or  three  dimensional  analysis  (See  Section  13.,  Part  IV): 

Next,  enter  a  2  If  a  two-dimensional  analysis  Is  desired. 

Enter  a  3  If  a  three-dimensional  analysis  Is  desired. 
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V.  Input  P|t> 


3)  Nowil  drag  fore as  (See  Section  1?.,  Port  IV):  Enter  a  1  to 
Include  these  forces.  Enter  e  0  to  neglect  these  forces. 

4)  Tangential  dreg  forces  (See  Section  17. ,  Pert  IV):  Enter  a  1 
to  Include  these  forces.  Enter  a  0  to  negleet  these  forces. 

5)  Added  Bass  force*  (See  Section  17. ,  Part  XV):  Enter  a  1  to 
Include  these  forees.  Enter  a  0  to  neglect  these  forees. 


A)  Bouyancy  forces  (See  Section  17.,  Part  XV):  Enter  a  1  to 
Include  these  forces.  Enter  a  0  to  negleet  these  forces. 

7)  Gravity  forces:  Enter  a  1  to  Include  these  forces.  Enter 
a  0  to  neglect  these  forces. 


8)  Vlecoscltles  and  fluid  aass  densities  (See  Section  17.,  Part  IV): 
Enter  a  0  to  use  the  default  values  coded  In  CABLE.  See  Table  I. 
If  values  different  froa  these  default  values  are  desired, 
enter  a  1.  Then  four  lines  of  data  Bust  be  entered,  providing 
the  following: 


I)  Vlscosclty  of  the  first  fluid  (usually  air)  Clb.sec./ft2 
orli.sec./B^'). 

II)  Bass  density  of  the  first  fluid  (slug/ft3  or  kg/n3). 

English  Metric 


Viscosity  (air) 


Bass 

Density 


(air) 


Viscosity  (eater) 


Bass 

Density 


(water) 


3.7188  x  10~7(lb-sec)/ft2 
2.378  x  1(T3  slug/ ft* 
3.516  x  10”5(lb-sec)/ft2 
1.9856  slug/ft3 


1.7802  x  10~5(N-sec)/m2 
1.2252  kg/a 3 
1.6831  x  10"3(N-sec)/n2 
1.0231  x  103  kg/»3 


TABLE  X.  Default  Values  for  Fluid  Vlscoseltles 
and  Bass  Densities 
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V«  Input  Data 


111) 

1v) 


Vlscosclty  of  the  second  fluid  (usually  water)  (lb.sec./ft^ 
or  N.sec/*2). 

Hass  danslty  of  tha  second  fluid  (slug/ft3  or  kg/a3). 


9)  Fluid  valocltlaai  CABLE  has  tha  option  of  allowing  tha  fluids 
to  Hava  constant  streaa  valocltlas  In  any  dlractlon.  .  To  use 
this  option  antar  a  1.  Next,  enter  two  Unas  of  data  providing: 
tha  X,  Y,  and  Z  coaponants  ralatlva  to  Inartla  spaca  of  tha 
valocltlas  of  tha  two  fluids  (fluid  1  (air)  on  tha  first  llna 
and  fluid  2  (water)  on  tha  sacond  llna). 

To  decllna  this  option,  antar  a  0.  This  naans  tha  fluids  will 
have  zero  velocity. 


10)  fluid  Interface  L.v.l:  CABLE  hs.  th.  option  of  .Hoping  th.  fluid 

Int.rfsc  Hvtl  Cth.t  Is,  th.  u.t.r  l.v.1)  to  b«  locpt.d  st  ,«  n„„n 
wrtlct  dlstanc.  b.lou  th.  mn  .hip  fr*».  To  u»  this  option  Mt.r 
.  1.  Then  on  th.  n.xt  tin.  of  Input,  Htir  th.  vortical  dlitancc 
Cposltlv.  dounu.rd)  fro.  th.  ..an  .hip  fru».  To  diclln*  this  option 
wit.r  •  0.  This  Mins  th.  fluid  Int.rfw.  »d  th.  man  ship  frut 
will  ba  on  tha  saae  laval. 

Tha  sat  of  options  data  will  consuae  batwaanlOand  17  lines,  depending 
upon  whathar  vlscoscltles,  densities,  -aloe It las,  and  the  fluid  Interface  level 
are  specified  by  tha  user. 


3.  Nuabar  of  Links 


Following  tha  Input  on  user  options,  antar  tha  nuaber  of  cable  links 
on  a  single  line. 
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V.  Input  Diti 


I 


ft 

4.  Individual  Link  Data 


The  next  input  data  la  a  aeries  of  data  seta  describing  the  physical 
and  geometrical  properties  of  each  link,  as  follows: 

1)  Link  number  (See  Section  6.,  Part  IV). 

2)  Label:  Any  name  up  to  20  characters  in  length. 

3)  Link  Reference  Point  Label  (or  "joint"  label)  (See  Section  4. * 

Part  IV):  Any  name  up  to  20  characters  in  length. 

4)  Link  number  of  adjacent  lower  numbered  link. 

5)  Link  weight  density  (weight  per  unit  length*  lb/ft  or  N/m). 

6)  Link  length  (ft  or  m). 

7)  Link  diameter  (in  or  cm). 

8)  Initial  configuration  of  the  link.  This  line  will  generally 
contain  three  numbers  representing  y,  a*  and  8  the  initial 
angles  (in  degrees)  of  the  link.  (See  Sections  11.  and  12. » 

Part  IV). 

If  the  two-dimensional  analysis  is  selected*  only  one  initial 
link  angle*  0,  is  entered. 

9)  Initial  motion  of  the  link.  This  line  will  generally  contain  three 
numbers  representing  y,  «*  and  0  the  derivatives  of  the  initial  angles 
of  the  link  (In  degrees/second).  If  the  two-dimensional  analysis  is 
selected*  only  one  initial  link  angle  derivative*  0*  is  entered. 

Nine  lines  of  data  are  required  for  each  link. 
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V.  Input  Data 


5.  Number  of  Towed  Spheres 

The  link  data  la  followed  by  analogoua  data  for  the  towed  sphere*. 

The  first  line  of  towed  sphere  data  contains  the  nuaber  of  towed  spheres. 

6.  Individual  Towed  Sphere  Data 

The  next  Input  data  Is  a  series  of  data  sets  describing  the  physical 
and  geometrical  properties  of  each  towed  sphere  as  follows: 

1)  Towed  sphere  nuaber  (See  Section  6.,  Part  IV). 

2)  Towed  sphere  label  Cup  to  20  characters). 

3)  Link  nuaber  of  the  adjacent  towing  link. 

4)  Towed  sphere  weight  (lb  or  N). 

5)  Towed  sphere  disaster  (In  or  ca). 

Five  lines  of  data  are  required  for  each  towed  sphere. 

7.  Potion  of  the  Wean  Ship  Frawe  and  the  System  Reference  Point 

The  aean  ship  frame  may  have  planar  motion  relative  to  the  Inertia 
frame  R.  Two  variables  representing  the  forward  (or  reverse)  speed  and 
the  turning  angle  are  used  to  describe  this  motion.  Also,  the  systea 
reference  point  d  may  have  arbitrary  motion  relative  to  the  mean  ship 
fraae,  requiring  an  additional  three  variables.  (See  Section  7.,  Part  IV.) 

For  the  Input  data  these  variables  are  Identified  by  the  following 
Integers: 
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V.  Input  Data 

>  • «  ■  * 

i  • 

1  -  Forward  displacement  (x)  of  reference  point  Q  relative  to  the 

Mean  ship  frame. 

2  -  Starboard  displacement  (y)  of  reference  point  Q  relative  to  the 

Mean  ship  frame. 

3  -  Vertical  (downward)  displacement  (z)  of  reference  point  Q 

relative  to  the  mean  ship  frame. 

4  -  Forward  (or  reverse*  If  negative)  speed  of  the  mean  ship  frame 

relative  to  the  Inertia  frame  R.  Note:  If  a  precoded  function 
1  fi  used  to  specify  this  motion,  the  mean  ship  frame  Is  assumed 
to  coincide  with  the  Inertial  frame  at  t  *  0.0. 

-4  -  _ .... 

5  -  Turning  angle  of  the  mean  ship  frame  relative  to  the  inertia 

frame  R.  (The  turning  angle  Is  positive  for  :  tarboard  turning.) 

The  data  describing  these  variables  may  be  specified  In  two  ways: 

1)  by  using  precoded  functions  and  2)  by  using  acceleration  profiles. 

(See  Sections  8.  and  9.,  Part  IV.)  The  precoded  functions  are  of  the 
form: 


f  *  f q  +  Ae^>tcos(pt  +  $)  (16) 

where  fg,  A,  b,  p*  and  $  are  constants.  The  acceleration  profiles  are 
described  In  Section  8.  of  Part  IV  of  the  Manual. 

Th*  specific  input  data  requirements  are  as  follows:  First,  enter  the 
number  of  these  variables  which  are  to  have  nonzero  values  during  the  motion. 
Then,  for  each  of  these  variables  enter  Its  Identifying  integer  (see  above) 
in  column  1.  In  column  2  enter  a  1  or  a  2  depending  on  whether  a  precoded 
function  or  an  acceleration  profile  is  desired. 
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V.  Input  Data 


If  a  precoded  function  It  selected,  enter  the  values  of  1q,  A,  b,  p, 
and  $  on  the  next  line. 

Alternatively,  If  an  acceleration  profile  Is  selected,  enter  the 
number  of  data  points  Cup  to  25)  on  the  next  line.  This  line.  In  turn, 
snjst  then  be  followed  by  the  same  number  of  data  lines  as  there  are  data 
points  (one  line  per  point).  The  first  of  these  lines  contains:  1)  the 
Initial  time,  11)  the  Initial  acceleration.  111)  the  Initial  speed,  and 
1v)  the  Initial  displacement.  The  subsequent  lines  contain  two  numbers: 

1)  the  time  and  11)  the  acceleration  value. 

If  several  variables  have  specified  motion,  their  data  Is  included  In 
serial  fashion. 

8*  Integration  Options 

The  next  data  is  a  set  of  four  numbers  setting  the  parameters  of  the 
Integration  subroutine  RKGS.  These  numbers  are  entered  on  a  single  line 
as  follows:  1)  Integration  starting  time,  11)  Integration  ending  time, 

111)  Integration  time  Increment,  and  tv)  Upper  bound  on  the  error. 

Note:  If  the  Integration  time  Increment  Is  Input  as  a  negative  value, 
say  -I,  then  CABLE  will  automatically  reset  the  Increment  to  (2)”1. 

9.  Output  Options 

1)  The  next  line  contains  the  time  increment  for  printing  the  computed 
data.  Since  printing  this  data  at  each  integration  step  can  produce 
an  excessive  amount  of  printed  output,  the  print  time  Increment  Is 
usually  greater  than  the  Integration  time  increment. 

2)  At  the  beginning  of  the  output  data?  CABLE  lists  a  copy  (or  "echo") 
of  the  Input  data.  The  extent  of  the  listing  of  computed  data  Is 

a  user  option.  To  exercise  this  option,  enter  an  Integer  (0,1, 2, 3, 4) 
In  column  1  according  to  the  following  printing  priority: 


V.  Input  Data 


0  -  Print#  all  computed  data:  Posit  Ion*  velocity*  and  accalaratlon 
of  tha  Man  ship  fraMt  Position*  velocity*  and  acceleration  of  Q 
the  systea  reference  polntt  Orientation  angles  and  their  deriva¬ 
tives  for  each  link  (In  degrees*  degrees/second,  and  degrees/ 
second/second) {  Position*  velocity*  and  acceleration  of  the  link 
and  towed  sphere  pass  centers  relative  to  Inertia  space  and  rela¬ 
tive  to  the  Man  ship  fraMt  Position  velocity,  and  acceleration 
of  the  connecting  Joints  relative  to  the  Inertia  space  and  relative 
to  the  Man  ship  fraMt  Cable  tension  at  the  systea  reference 
point  Q;  and  Cable  tension  at  the  connecting  joints* 

1  -  Prints  as  above  for  the  option  Integer  0*  except  for  the  cable 

tension  at  the  connecting  joints. 

2  -  Prints  as  above  for  the  option  Integer  1,  except  for  the  cable 

tension  at  the  systea  reference  point. 

3  -  Prints  as  above  for  the  option  Integer  2*  except  for  the  position 

velocity,  and  acceleration  of  the  joints. 

4  -  Prints  as  bove  for  the  option  Integer  3,  except  for  the  position, 

velocity,  and  acceleration  of  the  mass  centers. 
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VX.  EXAMPLE  INPUT  DATA 


For  an  Illustration  of  the  Input  data  consider  the  cable  link  systea 
shown  In  Figure  15. 


>  Figure  15.  Example  Cable  Systea 

Let  the  eleven  cable  links  be  Identical  having  length  2  ft.,  dlaaeter  0.25  In., 
and  weight  density  of  0.2  1b./ft.  Let  the  three  towed  spheres  have  dlaaeters 
1.0  In.,  1.5  In.  and  2.0  In.  with  weights  of  0.5  ob.,  1.688  lb.,  and  4.0  lb. 

t  Let  the  aean  ship  fraae  have  a  forward  aotlon  defined  by  the  acceleration 

i 

profile  shown  In  Figure  16.  Let  the  Initial  forward  ship  speed  be  5  ft/sec., 
and  let  the  Initial  dlsplaceaent  be  zero. 


Figure  16.  Acceleration  Profile  for  Forward  Mean  Ship  Frame  Motion 
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Let  the  systea  reference  point  (point  A)  have  oscillatory  Motion  relative 
to  the  aean  ship  fraae  given  by  the  expression 

z  ■  0.5  cos  *  t 

Let  the  fluid  level  be  3  ft.  below  the  aean  ship  fraae. 

Let  the  links*  spheres*  and  joints  be  labelled  as  shown  r,  Figure  15. 

Let  the  systea  be  aovlng  In  still  air  and  cala  water.  Let  gravity  and 
fluid  forces  be  acting  and  let  the  air  and  water  viscosities  and  densities  be 
the  saae  as  the  defeult  values  In  CABLE. 

Finally*  to  Illustrate  the  Input  of  Initial  conditions*  let  the  Initial 
configuration  of  the  cable  systea  be  In  the  X-Z*  plane.  Let  link  1  aake  an 
angle  of  45°  with  the  horizontal  and  let  link  2  aake  an  angle  of  30°  with  link 
1  as  shown  In  Figure  15.  Let  the  Initial  angles  of  all  other  links  all  be 
zero  (even  though  they  are  not  depicted  that  way  In  Figure  15.) 

The  Input  data  for  a  three-dimensional  analysis  of  the  notion  Is  then 
as  follows:  (Note:  The  coaaents  In  parentheses  are  for  Illustration  only 
and  should  not  aooear  In  the  noraal  FORTRAN  Input.) 

Line  Coluan  Nuaber 

1  2 

123456789012345678901234.  . . 

EXAMPLE  INPUT  DATA 

1 
3 
1 
1 
1 
1 
1 
0 
0 
1 

3.0 

1  2 

123456789012345678901234. . . 


(English  Units) 

(3-D  analysis) 

(Include  noraal  drag) 

(Include  tangential  drag) 
(Include  added  mass  forces) 
(Include  buoyancy  forces) 
(Include  gravity  forces) 

(use  default  fluid  properties) 
(fluids  are  at  rest) 

(specify  nonzero  fluid  Interface 
level) 
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1  2  3 

123456789012345678901234. ■ •  0 

11 
1 

LI 
A 

0 

0.2 
2.0 
0.25 

0.0  0.0  45.0 

0.0  0.0  0.0 

2 
L2 
B 

1 

0.2 
2.0 
0.25 

0.0  0.0  "*30.0 

0.0  0.0  0.0 

3 
L3 
C 

2 

0.2 
2.0 
0.25 

0.0  0.0  0.0 

0.0  0.0  0.0 

4 
L4 
D 

3 

0.2 
2.0 
0.25 

0.0  0.0  0.0 

0.0  0.0  0.0 

1  2 

123456789012345678901234. . . 


(nuaber  of  cable  links) 

(link  body  nuaber) 

(link  label) 

(reference  point  label) 

(lover  link) 

(wolght/ length) 

(length) 

(disaster) 

(Initial  orientation) 

(Initial  angle  derivatives) 

(next  link  body  nuaber  and  data) 


(next  link  body  nuaber  ard  data) 


(next  link  body  number  and  data) 


3 

0 
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1  2  3 

123456789012345678901234...  0 

5 

L5 

E 

4 

0.2 

2.0 

0.25 

0.0  0.0  0.0 

0.0  0.0  0.9 

7 

16 

C 

2 

0.2 

2.0 

0.25 

0.0  0.0  0.0 

0.0  0.0  0.0 

8 
L7 

F 

7 

0.2 

2.0 

0.25 

0.0  0.0  0.0 

0.0  0.0  0.0 

9 
L8 
G 

8 

0.2 

2.0 

0.25 

0.0  0.0  0.0 

0.0  0.0  0.0 

1  2  3 

123456789012345678901234...  0 


(next  link  body  nuabar  and  data) 


Cnaxt  link  body  nuabar  and  data) 


Cnaxt  link  body  nuabar  and  data) 


(naxt  link  body  nudbar  and  data) 


un*  coluan  lumber 

1  2  3  • 

123456789012345678901234...  0 

11 

0.2 

2.0 

0.25 

0.0  0.0  0.0 

0.0  0.0  0.0 

13  (next  link  body  nudber  and  data) 


L11 

I 

12 

0.2 

2.0 

0.25 

0.0  0.0  0.0 

0.0  0.0  0.0 

3 
6 

51 

5 

0.5 

1.0 

10 

52 

9 

1.688 

1.5 

14 

53 

13 

4.0 

2.0 

2 

31 


0.0 

0.5 

0.0 

3.1416 

42 

0.0 

0.0 

5.0 

0.0 

1  2  3 

123456789012345678901234. . . 0 


(nunber  of  towed  spheres) 

Cipher*  body  nuaber) 

(label) 

(body  nuaber  of  lower  link) 

(weight) 

Cdlaaeter) 

(next  spherical  body  nuaber  and  data) 


(next  spherical  body  nuaber  and  data) 


(nuaber  of  specified  variables) 
(variable  3  Is  specified  with  a 
precoded  function) 

(variable  4  Is  specified  with  an 
acceleration  profile) 
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3 

0 


Line  Column  NiMfctr 
1  2 

123456789012345678901234. . . 

5.0  2.0 

15.0  0.0 

0.0  20.0  0.25  0.001 

0.5 
0 

1  2  3 

123456789012345678901234...  0 


(Integration  parameters) 

(tlaa  Increment  for  printing  data) 
(output  option— print  all  data) 


APPENDIX 


SUMMARY  OF  INPUT  DATA  REQUIREMENTS 


This  Appendix  to  the  Manuel  provides  an  abbreviated  summary  of  the 
Input  data  requirements  for  a  quick  reference. 

1.  Heading  or  Title:  Heading  of- up  to  80  characters. 

2.  Options:  As  follows: 

1)  Units* 

Enter  1  for  English  units. 

Enter  2  for  metric  units. 

2)  Two  or  three  dimensional  analysis: 

Enter  2  for  two-dimensional  analysis. 

Enter  3  for  three-dimensional  analysis. 

3)  Normal  drag  forces: 

Enter  1  to  Include  tte  forces. 

Enter  0  to  neglect  the  forces. 

4)  Tangential  drag  forces: 

Enter  1  to  Include  the  forces. 

Enter  0  to  neglect  the  forces. 

5)  Added  mass  forces: 

Enter  1  to  Include  the  forces. 

Enter  0  to  neglect  the  forces. 

6)  Buoyancy  forces: 

Enter  1  to  Include  the  forces. 

Entar  0  to  neglect  the  forces. 
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’  SUMMARY  Of  INPUT  DATA  REQUMCTPTTS 


I 

i 

7)  Gravity  forces: 

Enter  1  to  include  the  forces. 

Enter  0  to  neglect  the  forces. 

8)  Viscoscities  and  fluid  mass  densities; 

Enter  0  to  use  the  default  values. 

Enter  1  to  use  substitute  values  entered  on  four 
subsequent  data  lines  as  follows: 

i)  Viscoscity  of  the  first  fluid 
ii)  Mass  density  of  the  first  fluid 
ill)  Viscoscity  of  the  second  fluid 
1v)  Mass  density  of  the  second  fluid 

9)  Fluid  velocities: 

Enter  0  for  aero  fluid  velocities. 

Enter  1  for  constant  fluid  velocities.  Then,  on 
‘  the  next  two  lines,  enter  for  the  two  fluids, 

the  X,  Y,  and  2  components  of  the  velocities. 

I 

10)  Fluid  Interface  Level: 

Enter  a  0  for  a  aero  interface  level. 

Enter  a  1  for  a  nonaero  interface  level. 

Then,  on  the  next  line,  enter  a  distance. 

3.  Number  of  links:  Enter  number. 

4.  Individual  link  data:  As  follows: 

1)  Link  number 

2)  Label:  Up  to  20  characters 

3)  Link  reference  point  label:  Up  to  20  characters 

4)  Link  number  of  adjacent  lower  numbered  link 

5)  Link  weight  density  (per  unit  length) 

6)  Link  length 

7)  Link  diameter 
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8)  Initial  configuration  of  the  link:  Enter  y,  a,  and  B 

Initial  angles  (In  degrees)  for  three-dimensional 
analysis.  Enter  0  for  two-dimensional  analysis. 

9)  Initial  notion  of  the  link:  Enter  y,o»  and  0  the  derivatives  of 
of  the  Initial  angles  (In  degrees/ second)  for  three-dimensional 
analysis.  Enter  0  for  two-dimensional  analysis. 

5.  Number  of  towed  spheres:  Enter  number. 

6.  Individual  towed  sphere  data:  As  follows: 

1)  Towed  sphere  number 

2)  Towed  sphere  label  (up  to  20  characters) 

3)  Link  number  of  adjacent  towing  link 

4)  Towed  sphere  weight 

5)  Towed  sphere  diameter 

7.  Notion  of  the  mean  ship  frame  and  the  system  reference  point: 

Enter  data  only  for  non-zero  variables  among  the  following: 

Identification  Number  Variable 

1  Forward  displacement  of  reference 

point  Q  relative  to  the  mean  ship  frame. 

2  Starboard  displacement  of  reference 

point  Q  relative  to  the  mean  ship  frame. 

3  Vertical  (downward)  displacement  of 

reference  point  0  relative  to  the  mean 
ship  frame. 

4  Forward  speed  of  the  mean  ship  frame. 

5  Turning  angle  (positive  starboard)  of 

the  mean  ship  frame. 

Enter  the  variable  Identification  number  In  column  1.  In  column  ?. 
enter:  1  for  preeoded  function*  ?  for  acceleration  profile. 
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Zf  a  preceded  function  It  talactad  values  of  fg,  A,  B,  p,  and  4 
on  the  subsequent  line*  (See  Equation  (16).) 

If  an  acceleration  profile  Is  selected,  enter  the  following  data  lines 

1)  Number  of  points  (up  to  25) 

2)  Initial  timet  Initial  acceleration!  Initial  speedt 
Initial  displacement  (all  on  one  line). 

3)  Time;  acceleration  (one  line  for  each  remaining  point). 


8-  integration  Options;  Enter  on  a  single  line:  1)  starting  timet 
11)  ending  timet  111)  time  Increment!  and  1v)  error  bound. 

9.  Output  Options: 

1)  Enter  the  printing  time  Increment 

2)  Enter  number  from  0  to  4  In  column  1  as  follows: 

0  -  Print  all  data. 

1  -  Print  all  but  tension  at  the  joints. 

2  -  Print  as  above  except  for  tension  at  the  reference  point  Q. 

3  -  Print  above  except  for  joint  kinematics. 

4  -  Print  as  above  except  for  mass  center  kinematics. 
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